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D , We present a modification of our previous family replicated gauge group model, 

' which now generates the Large Mixing Angle MSW solution rather than the ex- 

perimentally disfavoured Small Mixing Angle MSW solution to the solar neutrino 
^ \ oscillation problem. The model is based on each family of quarks and leptons having 

^ \ its own set of gauge fields, each containing a replica of the Standard Model gauge 

fields plus a (i? — L)-coupled gauge field. By a careful choice of the Higgs field 
gauge quantum numbers, we avoid our previous prediction that the solar neutrino 
mixing angle is equal order of magnitudewise to the Cabibbo angle, replacing it and 
the well-known Fritzsch relation with the relation Oc ~ (^o)"^''^ {f^d/ms)'^/^ . At 
the same time we retain a phenomenologically successful structure for the charged 
quark and lepton mass matrices. A fit of all the seventeen quark- lepton mass and 
mixing angle observables, using just six new Higgs field vacuum expectation values, 
agrees with the experimental data within the theoretically expected uncertainty of 
about 64%, i.e. it fits perfectly order of magnitudewise. 
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1 Introduction 



The first results on the charge current interactions from the Sudbury Neutrino Obser- 
vatory (SNO) collaboration [Q] have provided an important signal confirming the existence 
of the solar neutrino anomaly puzzle |||, |], ^, |^: SNO detected a flux of non-electron 
neutrinos, and Ut, among solar neutrinos after travelling from the core of the Sun to the 
Earth. Combination of the SNO results with previous measurements from other experi- 
ments reveals a confirmation of the standard solar model 0, whose predictions of the total 
flux of active neutrinos in the Sun agree with the SNO and Super-Kamiokande ||^ data. 
Furthermore, the measurement of the and hep solar neutrino fluxes shows no significant 
energy dependence of the electron neutrino survival probability in the Super-Kamiokande 
and SNO energy ranges. These results support the Large Mixing Angle MSW solution 
(LMA-MSW) rather than the Small Mixing Angle MSW solution (SMA-MSW) to the 
solar neutrino problem. 

Another important result on the solar neutrino problem, reported by the Super- 
Kamiokande collaboration p[, is that the day-night asymmetry data disfavour the SMA- 



MSW solution at the 95% C.L.. In fact, global analyses |T^, [n|, |T^, [T3| of solar neutrino 
data, including the first SNO results and the day-night effect, have confirmed that the 
LMA-MSW solution gives the best fit to the data and that the SMA-MSW solution is very 
strongly disfavoured and only accepted at the 3a level. The best fit values of the mass 
squared difference and mixing angle parameters in the two flavour LMA-MSW solution^ 
are Am| ^ 4.5 x 10"^ eV^ and tan^ 6q ^ 0.35. 

We have previously attempted to flt all the fermion - quark and lepton - masses 



and mixing angles including baryogenesis |T^, |T^] in a rather specific model without 



supersymmetry or grand unification. The model has the maximum number of gauge 
fields consistent with maintaining the irreduciblity of the usual Standard Model fermion 
representations, including three right-handed neutrinos. The predictions of this previous 
model are in order of magnitude agreement with all existing experimental data; however, 
only provided we use the SMA-MSW solution. But, for the reasons given above, the SMA- 
MSW solution is now disfavoured phenomenologically. So, in this article, we present a 
modified version of the previous model, which manages to accommodate the LMA-MSW 
solution for solar neutrino oscillations: all the fermion mass and mixing angle parameters 
are fitted within a factor of two, using 6 adjustable parameters. 

This article is organised as follows: in the next section, we define our notation for 
the charged fermion Yukawa coupling matrices, mass matrices and mixing angles. Then, 
in section 3 we review the family replicated gauge model. In section 4 we discuss the 
reasons for the modification of our model and the introduction of new Higgs fields. The 
calculation is described in section 5 and the results are presented in section 6. Finally, 
section 7 contains our conclusion. 



^The best fit parameter values for the LMA solution depend somewhat on the analysis method. 
However they do not change drastically from one two flavour analysis to another. We discuss the three 
flavour analyses in section |6.l|. 
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2 Charged fermion masses and their mixing angles 



In the Standard Model all fermions (apart from the neutrinos) get a mass via the 
electroweak spontaneous symmetry breaking - the Higgs mechanism. In extensions of 
the Standard Model containing right-handed neutrinos, the physical light neutrinos get a 
mass via the see-saw mechanism (see discussion of the see-saw mechanism in section |3.4D . 
The Higgs mechanism generates charged fermion mass terms from their Yukawa couplings 
in the Standard Model Lagrangian: 

~ -^chargcd-femiion-mass = Q l^u^wsUr + QlYo^wsDr + LiX^^^sEr + h.C. (1) 

Here ^ws is the Weinberg-Salam Higgs field, Ql denotes the three SU{2) doublets of 
left-handed quarks, Ur denotes the three singlets of right-handed up-type quarks and 
Yu is the three-by-three Yukawa coupling matrix for the up-type quarks. Similarly Yc 
and Ye are the Yukawa coupling matrices for the down-type quarks and charged leptons 
respectively. The SU (2) doublets ^ws and Ql can be represented as 2 component column 
vectors and we then define: 

^ws={ \ j)<fL (2) 

and 



= ( = D,) (3) 



L 



where [/^ are the CP conjugates of the three left-handed up-type quarks. After elec- 
troweak symmetry breaking the Weinberg-Salam Higgs field gets a vacuum expectation 
value (VEV) and we obtain the following mass terms in the Lagrangian: 

— '^^charged-fermion-mass = Ul Ur + Dr + Er + h.C. (4) 

where the mass matrices are related to the Yukawa coupling matrices and Weinberg-Salam 
Higgs VEV by: 

M = Y %^ (5) 
v2 

We have chosen the normalisation from the Fermi coupling constant so that: 

(0^,) = 246 GeV . (6) 

In order to obtain the masses from the mass matrices, M^/, and M^, we must 
diagonalise them to find their eigenvalues. In particular we can find unitary matrices, Vu 
for the up-type quarks, Vo for the down-type quarks and Ve for the charged leptons: 

V^M^MlVu = dis.g(mlmlml) (7) 

VlM^MlVr, = dmg(mlmlmt) (8) 

VIMeMIVe = diag (m2,mj,m2) (9) 

The quark mixing matrix is then defined with these unitary matrices as ||17||: 



Vo^M = V^V^ . (10) 



3 



3 Model with many quantum numbers 



We have already investigated a model |jTj, ^ which can predict not only quark and 



charged lepton quantities - masses and mixing angles - but also neutrino oscillations. 
This model has, as its back-bone, the property that there are generations (or families) 
not only for fermions but also for the gauge bosons, i.e. we have a generation (family) 
replicated gauge group namely 



X 

= 1,2,3 



{SMG, X f/(l),_,,J , (11) 



where SMG denotes the Standard Model gauge group = S'f/(3) x SU{2) x f/(l), x 
denotes the Cartesian product and i runs through the generations. For the prediction of 
the charged particle masses and mixings, the important part of the gauge group is the 
repetition of the Standard Model gauge group plus one extra U{1) called f/(l)/, where / 
denotes flavour [|18], |1^, ^ . But for the extension to neutrino masses and mixings using the 



see-saw picture, it is necessary to introduce a right-handed neutrino and a gauged B — L 
charge for each generation with the associated abelian gauge groups f/(l)s_i„i {i = 1, 2, 3). 
The just mentioned U{1) f abelian factor gets absorbed as a linear combination of the B—L 
charge and the weak hypercharge abelian gauge groups for the different generations. Note 
that this family replicated gauge group, eq. (0), is the maximal gauge group under the 
following assumptions: 



1 ) We only consider that part of the gauge group of Nature which acts non-trivially 
on the known 45 Weyl fermions of the Standard Model and the additional three 
heavy see-saw (right-handed) neutrinos. That is our gauge group is assumed to be 
a subgroup of ?7(48). 

2) We avoid any new gauge transformation that would transform a Weyl state from 
one irreducible representation of the Standard Model group into another irreducible 
representation: there is no gauge coupling unification. 

3) The gauge group does not contain any anomalies in the gauge symmetry - nei- 
ther gauge nor mixed anomalies. Note that otherwise the model becomes non- 
renormalisable. 



3.1 Gauge quantum numbers for the "proto" fermions at the 
fundamental scale 

In our model at the fundamental scale, which we take to be the Planck scale, there exist 
many bosons and fermions with practically all quantum numbers we can ask for. But 
most of the fermions have vector couplings, in the sense that they are described as Dirac 
particles from the Weyl point of view: they are combinations of left-handed and right- 
handed states with the same (gauge) quantum numbers. The left-over Weyl particles (in 
other word those without chiral partners) in our model are specified in more detail and 
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are actually assumed to form a system of three proto-generations, each consisting of the 
16 Weyl particles of a usual Standard Model generation plus one see-saw particle. In 
this way we can label these particle as proto-left-handed or proto-right-handed u-quark, 
(i-quark, electron etc. To get the quantum numbers under our model gauge group for 
a given proto-irreducible representation, we proceed in the following way: We note the 
generation number of the particle for which we want quantum numbers and we look up, 
in the Standard Model, what are the quantum numbers of the irreducible representation 
in question and what is the B — L quantum number. For instance, if we want to find the 
quantum numbers of the proto-right-handed strange quark, we note that the quantum 
numbers of the right-handed strange quark in the Standard Model are weak hypercharge 
y/2 = —1/3, singlet under SU{2) and triplet under SU{?>), while _B — L is equal to the 
baryon number = 1/3. Moreover, ignoring mixing angles, the generation is denoted as 
number i = 2. The latter fact means that all the quantum numbers for SMGi z = 1, 3 are 
trivial. Also the baryon number minus lepton number for the proto-generation number one 
and three are zero: only the quantum numbers associated with proto-generation two are 
non-trivial. Thus, in our model, the quantum numbers of the proto-right-handed strange 
quark are 1/2/2 = —1/3, singlet under SU{2)2, triplet under SU{3)2 and {B — L)2 = 1/3. 
For each proto-generation the following charge quantisation rule applies 

1 + 1 + 1 = (modi) , (12) 

where ti and di are the triality and duality for the z'th proto-generation gauge groups 
SU{3)i and SU{2)i respectively. 

Combining eq. (p!2D with the principle of taking the smallest possible representation of 
the groups SU{?>)i and SU{2)i, it is sufficient to specify the six Abelian quantum numbers 
yi/2 and {B — L)i in order to completely specify the gauge quantum numbers of the fields, 
i.e. of the Higgs fields and fermion fields. Using this rule we easily specify the fermion 
representations as in Table 1 (the representations of the Higgs fields will be given in 
subsection 5.1, where we present the fermion mass matrices). 

Note that each proto-generation gauge group SMGi x U {1)3-1,1 is a subgroup of 
S'O(IO), i.e. our gauge group eq. ( pJ]) is really a subgroup of 50(10)^. That means the 
i'th proto-generation has its own subgroup of S'O(10)j. However, we do not take the 
gauge fields of these SO{lQ)i to exist, except for those corresponding to the subgroups 

SMC, X U{l)B-^^. 



3.2 Breaking of the family replicated gauge group to the Stan- 
dard Model 

The gauge group ^=^2,3 {SMGi x U{1)B-L,i) is at first spontaneously broken down at one or 
two orders of magnitude below the Planck scale, by 5 different Higgs fields, to the gauge 
group SMG X U{1)b-l which is the diagonal subgroup of the original one: 

{ (f/, U,U)\U(Z SMG X U{1)b-l] C ,=^2.3 {SMGi x f/(l)s.^„) . (13) 
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Table 1: All ^7(1) quantum charges for the proto-fermions in the model. 
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We have to emphasize here that the gauge groups SMC and U{1)b-l act similarly on all 
three families, i. e. they are not any more family replicated gauge groups but correspond to 
the usual gauge group of the Standard Model and the usual baryon number minus lepton 
number. This diagonal subgroup is further broken down by yet two more Higgs fields 
— the Weinberg-Salam Higgs field ^ws and another Higgs field — to SU{3) x 

U{l)em- The vacuum expectation value (VEV) of the (pB-L Higgs field is taken to be 
about 10^^ GeV and is designed to break the gauged B — L quantum number. In other 
words the VEV {fpE-h) gives the see-saw scale. 

Let us stress that we have only one Weinberg-Salam Higgs field ^wsi i-^- it only has 
one irreducible representation in our family replicated gauge group. We freely use both 
^ws and its Hermitian conjugate ^ws, which means that we have no supersymmetry in 
the model preventing one or the other from giving masses to the quarks and leptons. 
Some of our predictions would be spoiled by introducing supersymmetry, because we 
need both a Higgs field and its Hermitian conjugate. With supersymmetry the number 
of Higgs fields and associated VEVs would have to be doubled; in the special case of the 
Weinberg-Salam Higgs field, this means introducing the unknown parameter tan (3. 
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3.2.1 Characterization of quark and charged lepton mass spectra 



An important prediction of our model depends on the strongly non-supersymmetric 
feature of there being only one Weinberg-Salam Higgs field, but it is independent of the 
details of the other Higgs fields which break our gauge group down to the Standard Model. 
This predicted feature is that corresponding diagonal matrix elements in each of the three 
charged mass matrices, Mu, M^, and even for the Dirac neutrino mass matrix, M^, 



are order of magnitudewise the same |118| , |T9[. The quantum number differences between 



the left- and right-handed Weyl fermions, between which a transition is needed to get 
these diagonal elements in our model, are yi/2 = ±1/2 and SU{2)i representation equal 
to doublet, the rest being trivial, where i is the proto-generation number of the diagonal 
element in question. Thus the quantum number violation needed, and therefore the order 
of magnitude resulting when all couplings are of order unity, will be the same for the 
diagonal element corresponding to a proto-family i in each of the four left-to-right mass 
matrices. My, Mj,, and (the Dirac neutrino mass matrix). 

The second and third family physical up-type quarks, t and c, get their masses from 
two off-diagonal elements in My which dominate the diagonal ones in our model 0, |l^ . 



So the above family degeneracy prediction then ends up becoming a prediction for the 
down-type quarks and leptons, simulating the simple SU{b) GUT prediction (m;, = m^, 
'^s = TT^ni t^d = t^e)i but we only get it with respect to order of magnitude. Thus 
our model can get the rough SU{5) mass predictions, without having to suffer from the 
problem of needing, say, an extra 45 Higgs at the Weinberg-Salam Higgs scale and thereby 
varying the Clebsch- Gordon coefficients so as to cope with, what is honestly speaking, 
sheer disagreement for the simplest SU{5) GUT. For the first family, in addition to the 
simulated GUT prediction, there is the degeneracy prediction that, when extrapolated to 
the Planck scale, m„ ^ rrie order of magnitude-wise. This is an example of a prediction of 
our model that is sensitive to it not being supersymmetric, because with supersymmety 
we would have two Weinberg-Salam Higgs fields and, with our philosophy that Higgs 
VEV's are likely to have their own order of magnitude, it would be difficult ever to get 
the prediction ^ mg. 

Another regularity predicted from our model is the "factorisation" of the quark mixing 
angles [1^, ^ 

Vub ^ Vus Vch . (14) 

This result mainly comes about because both Vuh and Vus contain, as a factor, similar Higgs 
field VEVs to take care of converting second family quantum numbers into first family 
ones. Really five of the eleven predictions of our model are made up from these general 
rules: four from the family degeneracy predictions and one from the above factorisation. 

3.3 Introduction of Right-handed Majorana neutrinos 

In order to explain the neutrino oscillations, we have introduced three very heavy right- 
handed neutrinos into our model, which are mass-protected by the Higgs field, (pB-L, at 
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an energy scale of about 10^^ GeV. We use the gauged B — L charge to mass-protect 
the right-handed neutrinos; in fact we use the total - diagonal - one because we break 
U{1)b~l,i X U{1)b-l,2 X U{1)b-l,3 3 U{1)b-l at a much higher energy scale, say about 
10^^ GeV. Another new Higgs field, x, "was also introduced in our previous see-saw model. 
This field plays the role of helping the VEV {(pB-i) to give non-zero effective mass terms 
for the see-saw neutrinos, by providing a transition between the right-handed tau-neutrino 
and the right-handed mu- neutrino. This transition coupling means that, with the new 
Higgs field we can obtain a large atmospheric neutrino mixing angle. 

However, unavoidably in the previous model, the solar mixing angle is in the region of 
the small mixing angle MSW (SMA-MSW) solution, i.e. the solar mixing angle and the 
Cabibbo angle are characterised by the same parameter, ^, of order 1/10. Furthermore 
the ratio of the solar neutrino mass squared difference to that for the atmospheric neutrino 
oscillations is given by fl^ without technical corrections ||2l| . On the other hand, with 



these technical corrections - "factorial factor corrections" - we could manage to make a 
mass squared difference ratio consistent with data for the SMA-MSW solutions : due to 
the presence of a Higgs field S, whose VEV is of order one in Planck units, there are many 
choices of the quantum numbers of the other Higgs fields that only change the number 
of occurrences of this field S in the fermion mass matrices. Moreover, one also has some 
freedom in the choice of the quantum numbers of the (ps-L field, which spontaneously 
breaks the gauged U{1)b-l group and thereby gives the see-saw scale (about 10^^ GeV). 
In this way, we managed to get the mass squared difference ratio to be of zeroth order in 
^ and rather to be given by 5*^/4, where the S field VEV is close to unity. However, the 
solar mixing angle could not be essentially changed, it remained of order ^, and thus only 
fits the SMA-MSW region. 



3.4 Neutrino masses and mixing angles 

The assumption of the existence of three right-handed Majorana neutrinos at a high 
scale0 gives rise to the addition of Majorana mass terms to the Lagrangian: 

- -^^neutrino-mass = i^L Ur + ^{WZY M L l^L + ^{T^Y M R Ur + k . C. 

\ijTZy MriL + h.c. (15) 



where 
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Here is the left-right transition mass term - Dirac neutrino mass term - and and 
Mr are the isosinglet Majorana mass terms of left-handed and right-handed neutrinos, 
respectively. 

^In the present model the right-handed neutrinos become massive by the action of the 4>b-l and x 
Higgs fields together with another new Higgs field p. These massive right-handed neutrinos would all 
have decayed and be washed out completely by the present epoch in the evolution of the Universe. 
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Due to mass-protection by the Standard Model gauge symmetry, the left-handed Ma- 
jorana mass terms, M^, are negligible in our model with a fundamental scale set by the 
Planck mass |2^. Then, naturally, the light neutrino mass matrix - effective left-left 
transition Majorana mass matrix - can be obtained via the see-saw mechanism p2| : 

M,s^M^M^\Mi;'f . (17) 

In the framework of the three active neutrino model, the flavour eigenstates Ua (a = 
e, I/, r) are related to the mass eigenstates z/j {i = 1,2,3) in the vacuum by a unitary 
matrix Kins, 

Wa) = Xl(^MNs)m Wi) ■ (18) 



Here Kins is the three-by-three Maki-Nakagawa-Sakata (MNS) mixing matrix |^ which 
is parameterised by 



K 



^ C13C12 

~C23Si2 - 
V S23S12 - Si3C23Ci2e 



Sl3S23Cl2e'"^-^ 



C13S12 
C23C12 ~ ■Si3S23Si2e 



-S23C12 - Si3C23Si2e' 



«di3 



Cl3'S23 
C13C23 



/ 



X 




(19) 





V 

where Ctj = cos 9ij and Sij = sin 9ij and ^13 is a CP-violating phase. Note that, due to the 
existence of Majorana neutrinos, we have two additional C7P-violating Majorana phases 
(p, ip, which are also included in the MNS unitary mixing matrix. 

In order to get predictions for the neutrino masses from the effective mass matrix, 
Meff, we have to diagonalise this matrix using a unitary matrix, Ves, to find the mass 
eigenvalues: 

KffMeffMjgKff = diag(m?, ml ml) . (20) 
With the charged lepton unitary matrix V^, eq. (^), we can then find the neutrino mixing 
matrix: 

KiNs = V;fr^^. • (21) 
Obviously, we should compare these theoretical predictions with experimentally measured 
quantities, therefore we define: 

(22) 
(23) 
(24) 
(25) 

Note that since we use the philosophy of order of magnitudewise predictions (see 
section ^ with complex order one coupling constants, our model is capable of making 
predictions for these three phases, the CP-violating phase ^13 and the two Majorana 
phases; put simply, we assume all these phases are of order 7r/2, i.e. essentially maximal 
CP violations. 



Ami 


= ml - 


- m\ 




III 

CO to 


- ml 


tan^ 6*0 


= tan^ 




tan^ ^atm 


= tan^ 


^23 
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3.5 No-Go theorem for large mixing angle in previous model 



We have traced the reluctance of our previous see-saw models to fit the large mixing angle 
MSW solution to the following feature of the Dirac neutrino mass matrix, M^: for every 
column, i.e. for all right-handed neutrinos in our notation, the first row elements - left 
electron ones - are smaller than the other matrix elements in the same column, by at 
least a factor of ^ ~ 1/10. With this property, we can indeed prove that the solar mixing 
angle cannot be bigger than of order ^, if we do not fine-tune the right-handed neutrino 
sector. 

Note that there is the possibility of getting a large solar neutrino mixing angle from the 
charged lepton sector, if it has big mixing relative to the proto-flavours [^. Our model, 
however, has an almost diagonal charged lepton mass matrix. Therefore we unavoidably 
obtain a small solar mixing angle, unless we re-arrange the Dirac neutrino mass matrix. 
That means that both the solar and atmospheric mixing angles must come from the Dirac 
neutrino sector in our present model. 

Our no-go theorem states that, provided there is essentially no mixing in the charged 
lepton sector and that the Dirac neutrino mass matrix, M^, obeys 



(M,^)i, < (M,^)2. e and (M, 



(M,^)3. e for ^ = 1,2, 3 



(26) 



the solar mixing angle cannot be larger than of order ^ ^ 1/10 in our previous model). 

This no-go theorem is even harder to circumvent if one has an 50(10) gauge group, 
because the up-type mass matrix is then very strongly related with the Dirac neutrino 
mass matrix and also the down-type mass matrix is similarly related with the charged 
lepton mass matrix |j25|. Really though it is necessary to exclude higher dimensional 
15*0(10) representations than say 10 for Higgs fields at the Weinberg-Salam Higgs scale. 



in order to obtain the identity of the mass matrices [pq] : 



M„ = M 



D 



Mo = 



(27) 



However, using this relationship, it is totally impossible to get a large solar mixing angle 
in the 50(10) model. 



4 Discussion of modification 



In order to get an LMA-MSW solution to the solar neutrino problem, we have to re-arrange 
the Dirac neutrino sector as we have already discussed in the previous section. We 
do this by the introduction of two new Higgs fields p and u, which replace 5* and ^. 

In our present model, we manage to make all the three elements of the first column 
(coupling to the first right-handed neutrino - the lightest one in our case) in the Dirac 
neutrino mass matrix roughly equal in order of magnitude, i.e. the (1, 1), (2, 1) and even 
(3, 1) matrix elements are made the same order of magnitude. The transition from the 
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second to the first column corresponds to a shift in the generation B—L quantum numbers, 
since it is given by the difference in the charges of the respective right-handed neutrinos. 
The new Higgs field, p, plays this role; more precisely the third power of p carries the 
quantum numbers required to make a transition form the first to the second column in 
the Dirac neutrino sector^]. 

In our new version of the model, the second and the third column in the Dirac neutrino 
mass matrix still obey the condition of our "no-go" eq. (p6D, which is actually very nice 
since it is needed to have a hope of getting a small CHOOZ angle ^13. If now all the 
elements in the first column would simply be obtained from the second by multiplication 
with {p^y as the quantum numbers at first suggest, this column would inherit the property 
of the first row element being small and we would not be able to get an LMA-MSW 
solution. However, we managed to get a need for the use of the p field in the matrix 
element (1,2) so that it has a factor of p^ in it, and then in the transition to the first 
column we get rid of the p^ factor rather than getting an extra factor of (p^)^. In this 
way we succeeded in making the ratio of matrix element (1, 1) to (2, 1) become bigger by 
a factor of than the ratio of (1, 2) to (2, 2). We really want the ratio of matrix element 
(1, 1) to (2, 1) to be of order unity, in order to obtain a large solar neutrino mixing angle. 
This is arranged by introducing the Higgs field u having a vacuum expectation value of 
about the same order of magnitude as p. 

The value of the Cabibbo angle corresponded to the VEV of ^ in the previous model. 
In the present model, it is given by the product up^ whose VEV should thus be of order 
of ^ ~ 1/10. From these considerations we can crudely estimate the VEVs of the new 
Higgs fields to be: a; ~ p ~ 1/3. 



5 Method of numerical computation 

A very important assumption in our model is that, at the Planck scale, we find a lot of 
different particles with many imaginable quantum numbers and having coupling constants 
which, when they are allowed, are complex numbers of order unity |2^ . This means that 
we assume essentially maximal CP violation in all sectors, including the neutrino sector. 
Since we do not know the exact values of all these couplings we are, in general, only 
able to make predictions order of magnitudewise. According to this philosophy |28| , 



we evaluate the product of mass-protecting Higgs VEVs required for each mass matrix 
element and provide it with a random complex number of order one as a factor. In this 
way, we simulate a long chain of fundamental Yukawa couplings and propagators making 
the transition corresponding to an effective Yukawa coupling in the Standard Model. In 
the numerical computation we then calculate the masses and mixing angles time after 
time, using different sets of random numbers and, in the end, we take the logarithmic 



^The quantum numbers of this Higgs field can be found in Table 
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average of the calculated quantities according to the following formula: 



{m)=exp[Y.^j . (28) 

Here (m) is what we take to be the prediction for one of the masses or mixing angles, rrii 
is the result of the calculation done with one set of random number combinations and N 
is the total number of random number combinations used. 

In order to find the best possible fit, we define a quantity which we call the goodness 
of fit (g.o.f.). Since our model can only make predictions order of magnitudewise, this 
quantity g.o.f. should only depend on the ratios of the fitted masses (mass squared differ- 
ences in the neutrino case) and mixing angles to the experimentally determined masses 
and mixing angles: 

Here (m) are the fitted masses and mixing angles defined in eq. (pSj) and mexp are the 
corresponding experimental values. The Yukawa coupling matrices are calculated at the 
fundamental scale, which we take to be the Planck scale. We use the first order renormal- 
isation group equations for the Standard Model to calculate the matrices at lower scales. 
Running masses are calculated in terms of the Yukawa couplings at 1 GeV (see section 
53). 



(29) 



5.1 Quantum numbers of the Higgs fields 

The model we present in this article has exactly the same gauge group and gauge quantum 
numbers for the fermions as in earlier versions 0, |15[ of our see-saw model. It is only 



the system of Higgs fields which have different gauge quantum numbers and they are 
presented in Table |[ The only essential change, even of the Higgs system, is that the 
fields u and p in the table replace the previous Higgs fields S and and take on different 
quantum numbers. As can be seen from Table |^, the fields u and p have only non-trivial 
quantum numbers with respect to the first and second families. This choice of quantum 
numbers makes it possible to express a fermion mass matrix element involving the first 
family in terms of the corresponding element involving the second family, by the inclusion 
of an appropriate product of powers of p and u. 

In previous versions of the model, this role of the p and u fields was played by the 
fields S and ^ with the quantum number combinations (ordered as in Table 2): 

S: (1.-1,0.-1,-1?) (30) 

(5.-6.0.0.3.-3) (31) 

It is with these quantum numbers that one gets the "no-go" situation for the LMA-MSW 
solution, since the solar neutrino mixing angle then satisfies 6q ^ ^ ^ Vus- It turns out 
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Table 2: All U{1) quantum charges of the Higgs fields. 



SMGi 



SMGo 
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1 

3 





that, fitting with these "old" quantum numbers, the vacuum expectation value of the 
field 5* is close to being unity in fundamental (Planck) units. Once the Higgs field S 
had a VEV of order unity, a large number of inessential modifications of the Higgs field 
quantum numbers became possible: one could add or subtract the quantum numbers of S 
to/from any of the other proposed Higgs fields a large number of times, without making 
any changes except in small details. Therefore, in the previous work, it was necessary to 
consider and make fits using these other possibilities. 

The new Higgs fields u and p turn out to have VEVs of the order of 1/3. So, in the 
present model, there are no fields with a VEV of the order of unity and thus no such 
ambiguities in the choice of Higgs field quantum numbers. In this way the "new" model 
escapes the "discrete" parameters of shuffling around the Higgs quantum numbers by 
multiples of those of S. So one now has a smaller amount of hidden fitting and a good 
fit should thus be considered a bit more impressive than in the previous model! The new 
model is in this way simplified compared to the old one. 

With the system of quantum numbers in Table ^ one can easily evaluate, for a given 
mass matrix element, the numbers of Higgs field VEVs of the different types needed 
to perform the transition between the corresponding left- and right-handed Weyl fields. 
The results of calculating the products of Higgs fields needed, and thereby the order of 
magnitudes of the mass matrix elements in our model, are presented in the following mass 
matrices: 



the up-type quarks: 
the down-type quarks 



V2 



V2 



(to^ypW^T^ W^T^ {W^yr I (32) 



u^pW{T^f W{T^f T3 I (33) 
\u^pW\T^Y W\T^Y WT 
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the charged leptons: 



the Dirac neutrinos: 




u\p^fW{T^f W{T^Y WT\ I (34) 



(pt)3|ytT2 w^T^ W^T\ (35) 

{p'^fW^T^X^ W^T^X^ W^T^ ] 

and the Majorana (right-handed) neutrinos: 

Mn ^ {<P^-J {pmxr ix^? X^ (36) 

V {p^fx^ x^ 1 ; 

In order to get the true model matrix elements, one must imagine that each matrix 
element is provided with an order of unity factor, which is unknown within our system 
of assumptions and which, as described above, is taken in our calculation as a complex 
random number, later to be logarithmically averaged over as in eq. ([281) . 

Note that the quantum numbers of our 6 Higgs fields are not totally independent. In 
fact there is a linear relation between the quantum numbers of the three Higgs fields W, 
T and x- 

Q^ = 3Qw-9Qt (37) 

where the 6 components of the charge vector Q correspond to the 6 columns of Table 
2. Thus the Higgs field combinations needed for a given transition are not unique, and 
the largest contribution has to be selected for each matrix element in the above mass 
matrices. 

Furthermore, there is another remark: the symmetric mass matrix - for the Majorana 
neutrinos - gives rise to the same off-diagonal term twice. The Feynman diagram for 
off-diagonal elements of the right-handed neutrino matrix is 




Thus to avoid overcounting we just have to multiply off-diagonal elements of the right- 
handed Majorana mass matrix by a factor of 1/2. However, in the Dirac mass matrix 
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columns and rows are related to completely different Weyl fields and, therefore, we do not 
need to worry about overcounting - the off-diagonal elements should not be multiplied by 
an extra factor of 1/2: 




The previous versions of our model predicted the Fritzsch relation 2P Vus = 0, 



md/mg (however only order of magnitude wise), provided that the VEV of the field S 
was of order unity, S* ~ 1. With the above mass matrices, this relation is now replaced 
by a relation involving the solar neutrino mixing angle: 



Vu. 



70 



m. 



(40) 



5.2 Renormalisation group running of coupling constants 

It should be kept in mind that the effective Yukawa couplings for the Weinberg-Salam 
Higgs field, which are given by the Higgs field factors in the above mass matrices multi- 
plied by some order unity factors (taken as random numbers), are the running (effective) 
Yukawa couplings at a scale very close to the Planck scale. Thus, in our calculations, 
we had to use the renormalisation group /3-functions to run these couplings down to the 
experimentally observable scale, i.e. ^ = 1 GeV where n is the renormalisation point. 
This is because we took the charged fermion masses to be compared to "measurements" 
at the conventional scale of 1 GeV. In other words, what we take as input quark masses 
are the current algebra masses, corresponding to running masses at 1 GeV, except for the 
top quark. We used the top quark pole mass instead: 

M, = ,„,(M) (l + li^^) , (41) 

where we set M = 180 GeV as an input, for simplicity. 

Using the notation in eq. (0), we can define the one-loop /3 functions for the gauge 
couplings and the charged fermion Yukawa matrices as follows: 



i« 2dgi 41 3 
dt 10 

^dg2 19 3 

lovr — — = fifn 

dt 16 
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-1^ 2dg3 - 3 

dt 

dY,r 3 _ ,A _ r_ /17 



levr^^ = \ {y,{y,)^ - r,(r,)t) + { - (i,? + + 8,|) } y, 
Ys = Tr(3yjr^ + 3rjr^ + rir^) , 

where t = In /i. 

In order to run the the renormahsation group equations down to 1 GeV, we use the 
following initial values: 



U{1) 
SU{2) 
SU{3) 



g,{Mz) = 0.462 , r7i(Mpi,nek) = 0.614 (43) 
g2iMz) = 0.651 , (72(Mpianck) = 0.504 (44) 
gsiMz) = 1.22 , (73(Mpianck) = 0.491 (45) 



5.3 The renormalisation group equations for the effective neu- 
trino mass matrix 

The effective light neutrino masses are given by an irrelevant, non-renormalisable 



dimension 5 term 31, 32 



ALeff = ^ C'' {eabHJD {e,M) , (46) 

where are left-handed Weyl lepton fields with the flavour index i and SU{2) weak 
isospin index a, and C^^ is a symmetric matrix of coefficients: 

C*J(^) = Y^^\M^y^Yl^ . (47) 

Here Mr is the Majorana mass matrix of the right-handed Majorana neutrinos, and Y^, 
is the Dirac neutrino Yukawa coupling matrix. 

The renormalisation group equations for the symmetric matrix, are given by 

167r2-— = (-3(72 + 2A + 2Ys) - - ijlf ij + C'^\ , (48) 



dt ^ ^' 2 

where A is the Weinberg-Salam Higgs self-coupling constant and 

Ys = Tx{'iYlY^ ^'iYlY^^YlY^) . (49) 

The mass of the Standard Model Higgs boson is given by M\ = X {4>ws)'^ and, for 
definiteness, we take Mh = 115 GeV thereby fixing the value of the Higgs self-coupling 
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A = 0.2185. These evolution equations can be rewritten using the elements of the light 
neutrino effective mass matrix Meg as running quantities: 

IGtt^^ = {-H + 2A + 2Ys) Meff - ^ [m,^ {Y.Ylf + (F^Fj) M,^) . (50) 

Note that the renormalisation group equations are used to evolve the effective neutrino 
mass matrix from the see-saw sale, set by {(pB-i) in our model, to 1 GeV. We should 
emphasize that we have used the approximation of ignoring the running of the Dirac 
neutrino Yukawa coupling constants between the Planck scale and the see-saw scale; 
however, this effect is small and so this approximation should be good enough for our 
order of magnitude calculations. 



6 Numerical results 

Using the three charged quark-lepton mass matrices and the effective neutrino mass matrix 
together with the renormalisation group equations, eqs. ( ^21) and (|50|), we made a fit to 
all the fermion quantities in Table 3 varying just 6 Higgs fields VEVs. We averaged 
over = 10,000 complex order unity random number combinations (see eq. |2^). These 
complex numbers are chosen to be the exponential of a number picked from a Gaussian 
distribution, with mean value zero and standard deviation one, multiplied by a random 
phase factor. We varied the 6 free parameters and found the best fit, corresponding to 
the lowest value for the quantity g.o.f. defined in eq. (pOj), with the following values for 
the VEVs: 

{(t)ws) = 246 GeV , = 1.64 x 10^^ GeV , {to) = 0.233 , 

(p) = 0.246 , {W) = 0.134 , (T) = 0.0758 , (x) = 0.0737 , (51) 

where, except for the Weinberg-Salam Higgs field and {(Pb-l), the VEVs are expressed in 
Planck units. Hereby we have considered that the Weinberg-Salam Higgs field VEV is 
already fixed by the Fermi constant. The results of the best fit, with the VEVs in eq. (|5lD, 
are shown in Table ^ and the fit has g.o.f. = 3.63. 

We have 11 = 17 — 6 degrees of freedom - predictions - leaving each of them with a 
logarithmic error of ^^3. 63/11 ~ 0.57, which is very close to the theoretically expected 
This means, in other words, that we can fit all quantities within a 
factor 1.78 ~ exp (^^3.63/11^ of the experimental value. However, we do not count the 
42 complex order unity random numbers entering the mass matrices eqs. (32)- (36) as 
parameters when discussing predictions, since we do not adjust them. We only use them 
as a calculational technique to avoid "unnatural" matrices which would be degenerate in 
our calculations, if we did not have random number coefficients. 

Unlike in older versions of the model, the first and second family sub-matrix of M^, is 
now dominantly diagonal. In previous versions of the model this submatrix satisfied the 
order of magnitude factorisation condition (Mo)i2 ■ {Mo)2i ~ {Md)ii ■ (Mo)22; thus the 



value 0.64 21 
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Table 3: Best fit to conventional experimental data. All masses are running masses at 
1 GeV except the top quark mass which is the pole mass. Note that we use the square 
roots of the neutrino data in this Table, as the fitted neutrino mass and mixing parameters 
(m), in our goodness of fit (g.o.f.) definition, eq. (|2^). 
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Vus 
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0.22 


Vcb 


0.027 


0.041 
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0.0025 
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9.5 X 10"^ eV^ 
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0.35 


tan 6a,tm 
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1.0 


tan^ 6*13 


4.8 X 10-2 


< 2.6 X 10-2 


g.o.f. 


3.63 





down quark mass rud received two contributions (off-diagonal as well as diagonal) of the 
same order of magnitude as the up quark mass m^. This extra off-diagonal contribution 
to md of course improved the goodness of the fit to the masses of the first family, since 
phenomenologically nid ~ 2 m^. However, in the present version of the model with the 
uj and p Higgs fields, the off-diagonal element {Mo)2i becomes smaller and we are left 
with a full order of magnitude degeneracy of the first family masses, even including the 
down quark. Furthermore, our expectation from section 4 that (u) ~ (p) ~ 1/3 tends 
to overestimate the first family masses. So the result of our fit, eq. (|51|), is to take (W) 
somewhat smaller than in our previous models and {u) ~ (p) < 1/3. Consequently our 
best fit values for the charm quark mass and the Cabibbo angle Vus are smaller than 
in our previous fits to the charged fermion masses |jl9|, |2ll]. Nonetheless, as mentioned 
above, our present best fit agrees with the experimental data within the theoretically 
expected uncertainty of about 64% and is, therefore, as good as can be expected from an 
order of magnitude fit. 

Experimental results on the values of neutrino mixing angles are often presented in 
terms of the function sin^ 29 rather than tan^ 9 (which, contrary to sin^ 29, does not have 
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a maximum at 9 = tt/A and thus still varies in this region). Transforming from tan^ ^ 
variables to sin^ 29 variables, our predictions for the neutrino mixing angles become: 



sin2 2^0 = 0.61 , (52) 
sin2 2^atm = 0.96 , (53) 
sin2 2^i3 = 0.17 . (54) 

We also give here our predicted hierarchical neutrino mass spectrum: 

mi = 4.9 X 10"^ eV , (55) 
m2 = 9.7 X 10^^ eV , (56) 
mg = 5.2 X 10"^ eV . (57) 

Compared to the experimental data these predictions are excellent: all of our order 
of magnitude neutrino predictions lie inside the 99% CL. border determined from phe- 
nomenological fits to the neutrino data, even including the CHOOZ upper bound. On 
the other hand, our prediction of the solar mass squared difference is about a factor of 
2 larger than the global fit data even though the prediction is inside of the LMA-MSW 
region, giving a contribution to our goodness of fit of g.o.f. ~ 0.14. Our CHOOZ angle 
also turns out to be about a factor of 2 larger than the experimental limit at 90% C.L., 
corresponding to another contribution of g.o.f. ~ 0.14. In summary our predictions for 
the neutrino sector agree extremely well with the data, giving a contribution of only 0.34 
to g.o.f. while the charged fermion sector contributes 3.29 to g.o.f.. 



6.1 CHOOZ angle and three flavour analysis 



The combination of the results from atmospheric neutrino experiments |^ and the 



CHOOZ reactor experiment constrains the first- and third-generation mixing angle to 
be small, i.e. the 3a upper bound is given by tan^ 6'i3 0.06. This limit was obtained from 
a three flavour neutrino analysis (in the five dimensional parameter space - 9q, ^13, 6'atm, 



Ami and Aml^^] 



using all the solar and atmospheric neutrino data and based on the 

2 A_2 ^ 



assumption that neutrino masses have a hierarchical structure, i.e. Arri'l <^ Aml^^^ 

However, the solar neutrino data in Table 3 come from a global two flavour analysis, 
which means that the first- and third-generation mixing angle is essentially put equal to 
zero, i. e. the dependence of ^13 on the solar neutrino parameters have been ignored. In 
principle we should, of course, fit to neutrino parameters from a three flavour analysis. 
Recently, global three flavour analyses were performed |3B, p^, and they showed a 



significant influence of the non-zero CHOOZ angle on the solar neutrino mass squared 
difference and mixing anglePI and vice versa: if the CHOOZ angle becomes far from zero 



''in iQ the relatively large solar neutrino mass squared difference lying in the LMA-MSW region 
(with the condition Attiq^IO"'* eV^), the solar mixing angle and the CHOOZ reactor experiment data 
were analysed using the three flavour analysis method. 
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then the solar mixing angle becomes smaller. This effect is more significant for the larger 
Am-Q values. Because of this correlation, our fit to the neutrino data is even somewhat 
better than that suggested by the g.o.f. value; even including the CHOOZ angle our 
neutrino fit is extremely good. 



6.2 CP violation 

We have fitted all the fermion masses and their mixing angles and therefore have predic- 
tions for the CKM and MNS mixing matrices, in the quark sector and in the lepton sector 
respectively, including CP violating phases of order unity. In this subsection, we will 
first consider the size oi CP violation in the quark sector and then the electron "effective 
Majorana mass" responsible for neutrinoless double beta decay. 

The Jarlskog invariant Jcp provides a measure of the amount oi CP violation in the 
quark sector |^ and, in the approximation of setting cosines of mixing angles to unity. 



is just twice the area of the unitarity triangle: 

Jcp = Vus Vcb Vub sin 6 , (58) 

where 6 is the CP violation phase in the CKM matrix. In our model the quark mass 
matrix elements have random phases, so we expect 6 (and also the three angles a, (3 
and 7 of the unitarity triangle) to be of order unity and, taking an average value of 
I sin (5 1 ~ 1/2, the area of the unitarity triangle becomes 

JcP~^K.KfeK6 . (59) 

Using the best fit values for the CKM elements from Table |^, we predict Jcp ~ 3.1 x 10~^ to 
be compared with the experimental value (2 — 3.5) x 10~^. Since our result for the Jarlskog 
invariant is the product of four quantities, we do not expect the usual ±64% logarithmic 
uncertainty but rather ±-\/4-64% = 128% logarithmic uncertainty. This means our result 
deviates from the experimental value by log |^Y^^^Tr6/1.28 = 1.7 "standard deviations". 

Another prediction, which can also be made from this model, is the electron "effective 
Majorana mass" - the parameter in neutrinoless beta decay - defined by: 



m, 



(60) 



where rrii are the masses of the neutrinos Vi and Uei are the MNS mixing matrix elements 
for the electron flavour to the mass eigenstates i. We can substitute values for the neutrino 
masses from eqs. ( p5| - p7| ) and for the fitted neutrino mixing angles from Table 3 into 
the left hand side of eq. (0). As already mentioned, the CP violating phases in the MNS 
mixing matrix are essentially random in our model. So we combine the three terms in 
eq. (^) by taking the square root of the sum of the modulus squared of each term, which 
gives our prediction: 

|(m)| ^ 3.1 X 10^^ eV . (61) 
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Although the Jarlskog invariant and the effective Majorana electron neutrino mass 
have been calculated from the best fit parameters in Table 3, it is also possible to calculate 
them directly while making the fit. So we have calculated Jcp and for = 10, 000 
complex order unity random number combinations. Then we took the logarithmic average 
of these 10, 000 samples of Jcp and \ {m)\ and obtained the following results: 

Jcp = 3.1 X 10"^ , (62) 
\{m)\ = 4.4 X 10"^ eV . (63) 

in good agreement with the values given above. 

We should mention here that our effective Majorana mass parameter eq. (|63|) , of 
course, respects the upper limit presented in ref. EDI . 



7 Conclusion 

We have developed an older version of our model, with the purpose of making it fit the 
experimentally favored LMA-MSW solution rather than the SMA-MSW solution for solar 
neutrino oscillations. In the older version, the magnitudes of the solar mixing angle ^0 and 
the Cabibbo angle 6c are both characterised by the VEV of the Higgs field ^ ~ 1/10 and 
thus the previous model could only be made compatible with the SMA-MSW solution. 
The required modification of the model was achieved by replacing the fields S and ^ in the 
previous model by another pair of Higgs fields: uj and p having non-trivial and opposite 
quantum numbers with respect to the family one and family two gauge groups, while 
having trivial family three gauge quantum numbers. In this way an excellent fit to the 
LMA-MSW solution is obtained. The price paid for the greatly improved neutrino mass 
matrix fit - the neutrino parameters now contribute only very little to the g.o.f. - is a slight 
deterioration in the fit to the charged fermion mass matrices. In particular the predicted 
values of the quark masses ma and rric and the Cabibbo angle Vus are reduced compared 
to our previous fits. However the overall fit agrees with the seventeen measured quark- 
lepton mass and mixing angle parameters in Table 3 within the theoretically expected 
uncertainty ||^ of about 64%; it is a perfect fit order of magnitude wise. 



It should be remarked that our model provides an order of magnitude fit /understanding 
of all the effective Yukawa couplings of the Standard Model and the neutrino oscillation 
parameters in terms of only 6 parameters - the Higgs field vacuum expectation values. So 
we can say that we fit all the parameters of the Standard Model and neutrino oscillations, 
except for the gauge coupling constants and the Higgs mass and its self- coupling. Actually 
we should note here that even the gauge coupling constants may be derived from order 
one quantities in the following sense: We postulate that the gauge couplings - and here, 
perhaps a bit arbitrarily at first, let us say that we take these to be gi, g2 and - are 
of order unity at the Planck scale. That means that the corresponding inverse a's (the 
inverse fine structure constants) are of order An = 12.5. Now, according to our model, the 
Standard Model gauge groups which we see experimentally are the diagonal subgroups of 
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a cross product of three replicas of the same mathematical group; one cross product factor 
for each family. The formula for calculating the inverse a for the diagonal subgroup, ol^s^ 
is: 



3 



a 



Ds = E . (64) 

i=l 

Thus the order unity assumption leads in our model to a Standard Model inverse a at 
the Planck scale being of the order of 3 ■ 12.5 = 37.5. This agrees very well with the 
experimentally measured Standard Model couplings, when extrapolated to the Planck 
scale using the renormalisation group equations of section 5.2, being respectively: ^ = 
55.5, = 49, a^^ = 54. This observation can be further elaborated and converted 
into an exact prediction by using the so-called Multiple Point Principle (MPP). The 



phase transition couplings used in this MPP can be taken as the definition of couplings 
being of order unity and, with this choice, it turns out that it is indeed the gi, g2 and 
that are really of order unity rather than, say, the ai, a2 and 0:3. 



In a similar way the application of the MPP, requiring degenerate minima in the 
Weinberg-Salam Higgs effective potential, can lead to a Higgs field self-coupling which is 
of a similar order to the squares of the giS, i.e. A ~ gf. 

Thus, including this sort of argument, our model gives an order of magnitude under- 
standing of all the Standard Model gauge, Higgs and Yukawa couplings, and even also of 
the beyond the Standard Model neutrino oscillation parameters. 
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